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Abstract 

|We examine the statistics of active scalar fluctuations in high-Rayleigh number fingering convection with high-resolution three- 
dimensional numerical experiments. The one-point distribution of buoyancy fluctuations is found to present significantly non- 
Gaussian tails. A modified theory based on an original approach by Yakhot (1989) is used to model the active scalar distributions as 
a function of the conditional expectation values of scalar dissipation and fluxes in the flow. Simple models for these two quantities 
^ Jiighlight the role of blob-Uke coherent structures for scalar statistics in fingering convection. 
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^ 1. Introduction 



Fingering convection is a peculiar convective flow, character- 
ized by a counter-gradient density transport, which is of inter- 
est for a wide range of fields, including stellar physics, metal- 
lurgy and volcanology(author?) It raises particular interest 
^ in oceanography (author?) l'^, '5], where warm, salty, and light 
. ^ waters floating above fresher, colder and denser waters create 
^^finger-favorable conditions. This situation is commonplace in 
(— I the thermocline of subtropical oceans, where finger-generated 
C!L(,diapycnal mixing is believed to affect the meridional overturn- 
ing circulation, and the uptake of heat and carbon dioxide (au- 
jthor?) 101 ■ Fingering convection occurs when two buoyancy- 
J> changing scalars with different diff'usivities are stratified in such 
OO a way that the least-diffusing one, if taken alone, would produce 
an upward, unstable density gradient, but the most-diffusing 
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pne reverses this tendency and produces a net downward den- 



• sity gradient. This set-up allows for a doubly-diffusive insta- 
bility, where infinitesimal perturbations to the initial stratifica- 
tion may undergo exponential growth. If a fluid parcel is dis- 
placed downward it looses the stabilizing, most-diffusing scalar 
r* at a faster rate than the destabilizing, least-diffusing one, be- 
. ^ pause of the difference in diffusivities. This results in an in- 
crease of the density of the parcel, which sinks at a lower depth, 
H where it loses even more stabilizing scalar A symmetrical ar- 2. The Simulations 



regime, the elongated, finger-like columns of the linear theory 
break down into shorter plumes, or blobs, as we shall call them, 
of rising and sinking fluid, having an aspect ratio much closer to 
one. Vigorous convection arises, dominated by the complicated 
dynamics of the interacting blobs (author?) In some 

instances a secondary instability disrupts the linear profiles of 
horizontally averaged temperature and saUnity, leading to the 
formation of the so-called "staircases" where finger zones are 
alternated with slabs of vertically well-mixed fluid (author?) 

Hi. 

While ample literature has been devoted to convection prob- 
lems with a single scalar, the Rayleigh-Benard set-up being the 
leading example, convection with two active scalars, such as 
fingering convection, has yet to be explored to the same ex- 
tent. In particular, knowledge on the statistical properties of 
scalar fluctuations can be useful for constraining and testing 
models of the convective scalar fluxes. In this paper we use 
high-resolution three-dimensional numerical simulations to ex- 
plore the high-Rayleigh number regime of fingering convection 
and we analyze and interpret the scalar fluctuation distributions 
in these flows using a novel formulation of a classical theoreti- 
cal approach. 



gument holds for a fluid parcel displaced upward. While the 
two scalars are transported along their gradients, a net counter- 
gradient buoyancy flux results from this mechanism. In the lin- 
ear theory maximum growth rate is attained at horizontal wave- 
lengths sufficiently small so that the difference in diffusivity of 
the two scalars is important, but not so small that the damping 
effect of viscosity becomes dominant (a few centimeters in the 
oceans), and at a vertical wavelength corresponding to the ver- 
tical domain extension (author?) |l5i|6[]. In the fully nonlinear 
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2.1. The equations of fingering convection 

As customary in this problem, we denote the least-diffusing 
scalar as salinity, S , and the most-diffusing one as temperature, 
T, even if the actual physical nature of the two scalars may 
be different. We consider a fluid layer of thickness d confined 
above and below by perfectly conducting, parallel, plane plates 
mantained at constant temperature and salinity. We bring the 
problem to a non-dimensional form by scaling temperatures and 
salinity with their plate differences Ar and ts.S , scaling lengths 
with the layer thickness, d, and using the haline diffusive time, 
Ts - d^/Ks, as a timescale, with ks the haline diffusivity. The 
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control parameters of the problem are the Prandtl, Lewis, ther- 
mal and haline Rayleigh numbers, defined as 

y KT gaATd' gpAS 
Pr - — , Le - — , Rt , "s = , 

respectively. In these expressions v is the kinematic viscosity, 
kt is the thermal diffusivity, g is the modulus of the gravity ac- 
celeration and a and fi are the thermal and haline linear expan- 
sion coefficients. From these the density ratio Rp - LeRj/Rs 
can be defined. Linear stability analysis shows that a necessary 
condition for fingering instability is 1 < Rp < Le. The density 
ratio controls directly the intensity of convection in the flow. 
Values close to one lead to fast instability and violent convec- 
tion. The haline diffusive time is the longest time scale in this 
problem. Far from marginality, it is often more convenient to 
use the convective time scale = {Pr Le Rs Y^^^ts ■ 

In the Boussinesq approximation, the resulting non- 
dimensional equations for fingering convection are: 




— -H U Vu 

dt 
dT 



-Vp + Pr Le \RsBi. + V^u] , 
LeV^T, 



dS_ 
Yt 



Vu 



0, 
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(2) 

(3) 
(4) 



where u - (u, v, w) is the solenoidal velocity field of the fluid, 
p is the pressure, and we have defined the buoyancy field 
B = RpT - S. The latter is the dynamically important lin- 
ear combination of the T, S scalar fields as it appears in the 
forcing term of the momentum equation ([1). Notice that the 
equations could be rewritten entirely in terms of buoyancy to- 
gether with some other linear combination of temperature and 
salinity. In this case buoyancy would play the role of the active, 
momentum-generating scalar, and the other scalar could be con- 
sidered almost passive, interacting only with buoyancy through 
a diffusive term. For example in the oceanographic literature 
a quantity called "spice" is sometimes defined (author?) ifioll 
as the linear combination of temperature and salinity which is 
maximally independent of buoyancy. 

We integrate the equations numerically with a code which is 
pseudospectral in the horizontal directions and with finite dif- 
ferences in the vertical, with a non-homogeneous vertical grid 
in order to better resolve the boundary layers and with a third- 
order fractional step method for time advancement (author?) 
lfl4l [Tsl [l6ll . Free slip boundary conditions are used at the top 
and at the bottom and a laterally periodic domain is assumed. 

While the density ratio represents the main control param- 
eter in this type of flows, in the sense that small variations 
of Rp determine large variations in the fluxes, we study the 
changes in the flow determined by changes in the vertical ex- 
tent of the domain, as encoded in the Rayleigh numbers. It 
has been suggested that simulations in a vertically periodic do- 
main should yield results comparable with those that could be 
observed in a physical situation where the vertical extension is 
many orders of magnitude larger than the scale of the fingers 






Figure 1: Rendering of buoyancy fluctuation iso-surfaces at tlie level 2.5(Tb, 
once statistical stationarity has been reached, in the center of the domain be- 
tween z = -0.2 and z = 0.2. The lateral panels show sections through the 
flow with a colormap saturating at Scrg. Yellow(blue) corresponds to posi- 
tive( negative) buoyancy fluctuations. 



(e.g. (author?) iQH). However, recent work (author?) Ill 

casts shadows on the meaningfulness of investigations in such 
a domain. This requires treating the magnitude of the Rayleigh 
numbers, and not only their ratio, as an independent free con- 
trol parameter of the problem. Accordingly here we use a con- 
figuration confined between upper and lower rigid plates, and 
push up the Rayleigh numbers while maintaining a fixed den- 
sity ratio Rp - 1.2, which is sufficient to guarantee vigorous 
convection. We perform a series of numerical experiments, fix- 
ing Pr - \0,Le - 3), and exploring the range Rs - 10** - 10". 
The vertical domain extension is -0.5 < z < 0.5. The ver- 
tical resolution ranges between - 145 (at Rg - 10**) and 
A^j. = 809 layers (at Rs - 10"). The horizontal resolution is 
maintained at = 192 grid points, while the lateral domain 
size is scaled in order to be approximately proportional to the 
horizontal linear instability scale (author?) |i6|], according to 
L = {Rsl\Q'^)-^^\ 

2.2. Statistics of scalar fluctuations 

Fig. [T] shows an isosurface rendering of the buoyancy fluc- 
tuation field for the simulation at - 10", in the center of 
the domain, after the solution has reached statistical stationar- 
ity. The flow is characterized by the presence of well defined 
buoyancy structures which transport a large fraction of the ver- 
tical buoyancy fluxes and lead to the vertical homogenization 
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Figure 2: Horizontally and time averaged scalar profiles for the simulation at 
Rg = lO" ; (in the inset) profiles of the standard deviation of scalar fluctuations 
at fixed z- The time averages extend over one hundred convective time units. 



of the statistical properties of the flo'wQ. Fig. |2] shows ver- 
tical profiles of horizontal averages of buoyancy, temperature 
and salinity, and of the variance of fluctuations respect to these 
profiles. Note the inverse boundary layers of buoyancy, which 
testify of a counter-gradient advective transport. An ample cen- 
tral region characterized by a uniform linear background and 
approximately z-independent variance is evident. This is the 
region of interest for many geophysical applications of finger- 
ing convection, including oceanographic problems, all charac- 
terized by high Rayleigh numbers and reduced influence of the 
vertical boundary conditions. Increasing the Rayleigh numbers 
in our simulations leads to a growth of the height of this central 
region. The presence of a vertical average background gradient 
of the scalars is characteristic of fingering convection, and it 
is consistent with laboratory experiments and observations (e.g 
(author?) Hi]) where approximately uniform vertical mean 
gradients are reported within the fingers zones. 

In the central region buoyancy fluctuations (respect to hor- 
izontally averaged buoyancy), vertical velocities, and convec- 
tive fluxes all vary according to approximate power laws with 
the Rayleigh number, as summarized in table [T] The table also 
reports two independent estimates of the characteristic scales 
of the structures in the buoyancy field: is the position of the 
first zero of the spatial autocorrelation function of buoyancy, 
computed keeping constant the coordinates y and z; Z,, is the cu- 
bic root of the average volume of the connected regions with 
\B'\ > Ictb'. Both estimates give similar results. 

Figs. Oa-c) compare the one-point probability distributions 
of temperature, salinity and buoyancy fluctuations, for haline 
Rayleigh numbers Rs = 10^ 10'', 10"', 10", computed at the 
center of the domain, -0.2 < z < 0.2, once statistical station- 
arity has been reached. Temperature and salinity present sub- 



' For example in the data in Fig. [T] points with buoyancy fluctuation larger 
than 2 standard deviations in modulus cover only 5% of the volume but carry 
37% of the total vertical buoyancy flux. 
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Figure 3: {a, b, c) Probability density functions of the temperature, salinity and 
buoyancy fluctuations respectively, normalized to unitary variance. Data for 
four different simulations at Rg = 10**, lO', lO'", lO" are reported, gathered 
over a time interval of one hundred convective times, once stationarity has been 
reached. For clarity the latter three curves are multiplied by two, four and 
eight units, respectively. Superimposed on the numerical data are expression 
(7) and its analogues for temperature and salinity, evaluated using the fits (|D 
(for temperature and salinity), ^9) (for buoyancy) and E(T\X)IX = X. Gaus- 
sian disfi'ibutions with unitary variance are reported for reference, using dashed 
lines. 
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Exponent: 
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0.39 


0.25 
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Table 1 : Statistics of buoyancy fluctuations as a function of haline Rayleigh 
number, computed at the center of the domain (-0.2 < z < 0.2), once sta- 
tistical stationarity has been reached. The standard deviations of buoyancy 
fluctuations and of vertical velocities, erg/ and crw, and the vertical convec- 
tive fluxes, {WB'), are reported; d and are two independent estimates of the 
size of the structures in the buoyancy field (see text for details). The rightmost 
column reports the Reynolds number of the buoyancy structures estimated as 
Re = (T\f/lxl(LePr). The bottom row reports exponents obtained fitting a power 
law scaling respect to Rs . 



Gaussian distributions at low Rs , which become roughly Gaus- 
sian at intermediate Rs and develop tails longer than Gaussian 
only at the highest Rs explored in this work. The distribution of 
buoyancy, instead, shows non-Gaussian tails already at the low- 
est Rs considered, which become very similar to exponentials 
at high Rs . 

The emergence of exponential tails in the active scalar dis- 
tribution is reminiscent of the classical transition to so-called 
"hard turbulence" for thermal convection, first reported in (au- 
thor?) OlSll . In that case buoyancy is determined by only one 
scalar (temperature), which develops well defined exponential 
tails at high Rayleigh number The behaviour of a passive scalar 
in homogeneous turbulence has been discussed in e.g. in (au- 
thor?) 1 1 1, 12] and also in that case non-Gaussian behavior is 
found. The present case is only partially similar to those just 
mentioned: no scalar is completely passive in fingering convec- 
tion, and the buoyant structures have a very different origin and 
morphology compared to single-scalar convection (finger-like 
blobs vs. turbulent plumes). 



3. A model for the fluctuation distributions 

3.1. Yakhot's approach 

The scalar probability distributions described above can be 
interpreted by adapting to the present case with two active 
scalars the theory developed by Yakhot(author?) 1. 1 3 .1 for the 
case of high-Rayleigh number thermal convection. An impor- 
tant pre-requisite of that theory is a constant, non-zero vertical 
gradient of the horizontally averaged density field, a situation 
which does not occur in the bulk flow of Rayleigh-Benard con- 
vection, but whose presence is, as we have discussed above, the 
hallmark of fingering convection. 

In order to proceed along these lines we express the buoy- 
ancy field as B(x,y,z, t) - B'(x,y,z, t) + GbZ, where Gb is the 
vertical gradient of the horizontally averaged buoyancy, and B' 
is the buoyancy fluctuation. The quantities T', S', Gj or Gs can 
be defined analogously. Combining Eq. (|2]i and (O we obtain 
an evolution equation for the buoyancy fluctuation, written in 
terms of buoyancy and temperature: 



dB' 
dt 



We multiply this expression by B'2"-i and denote with (■) a time 
and volume average. The volume averages are carried over a 
central portion of the domain, where the assumption of con- 
stant gradients holds (represented by the dashed lines of Fig. |2] 
for our simulations). Assuming a statistically steady state we 
can integrate by parts; the lateral boundary terms are zero for 
periodic (or no flux) boundary conditions; the top and bottom 
terms cancel each other because statitical stationarity implies 
that the vertical flux of all moments of B' must be the same at 
all heights. Dividing the resulting equation for a generic n by 
the equation for « = 1 we obtain 



(6) 



where we have defined the normalized buoyancy fluctuation, 

Xb - B' I (b'^) ^ , the normalized buoyancy flux, "Fb - 
wB' I {wB') and the normaUzed buoyancy dissipation rate, = 



«p(Lc-l)V7''+VB'] 



We refer to the latter quantity as 'buoy- 



fip(Lc-l)Vr'+VB'])' 

ancy dissipation', even if it may be negative somewhere in the 
fluid. In fact, the signature of doubly-diffusive convection is the 
ability of the second-order derivative terms in the equations to 
behave as sources of buoyancy fluctuations, rather than solely 
as sinks. But the large-scale, average effect of these terms re- 
mains that of sinks of variance: for n - \ Eq. (|6]l is the balance 
between the advective rate of extraction of buoyancy variance 
from the vertical gradients, and its diffusive dissipation rate at 
small scales. Note that this process cannot be described solely 
in terms of buoyancy: temperature or salinity fluctuations must 
appear explicitly in the buoyancy variance dissipation rate. 

In the following we interpret the quantities Xb, Tb, and;^'^ as 
stochastic variables and we make our last assumption, namely 
the equivalence of space-time and ensemble averages. We ob- 
serve that the probability density of Xb must have a compact 
support, because Eqs. (|2]l and (|3]l satisfy a maximum principle, 
and this forbids arbitrarily large buoyancy fluctuations; there- 
fore the averages that appear in Eq. (|6]l exist for any integer 
n>0. 

From Eq. it is possible to follow Yakhot's analy- 

sis(author?) llisll slavishly, obtaining an explicit expression for 
the probability density of Xb'. 



P(Xb) 



E(xb\O)P(0 ) 
E(xb\Xb) 



■ exp 



Jo 



xE(j(b\x) 



dx 



(7) 



+ U-WB' - wGb = Rp {Le - 1) V^r' + V^B'. 



(5) 



where £'( | ) denotes a conditional average. This is an exact 
relationship for the one-point probability density of buoyancy 
fluctuations which depends only on two unknown functions, 
namely the expected values of the buoyancy fluxes, E{'Fb\Xb), 
and of the dissipation, E(xb\Xb), conditioned on the buoyancy 
fluctuations. The constant f (0) is fixed by the normalization re- 
quirement of the density function. Far from the physical bound- 
aries it is reasonable to assume that the actual solution of the 
Boussinesq equations ([T]|4| will have the same up-down sim- 
metry of the equations themselves: It follows that P{Xb), and 
the unknown conditional expectations E(xb\Xb) and E{'Fb\Xb), 
are even functions of Xb- Furthermore, E{TBWi - 0, because 
Tb^Q where Xb = 0. 



Theoretical expressions for the distributions of normalized 
temperature and salinity fluctuations, or of their linear com- 
binations, such as spice, can be obtained following the same 
approach. This yields expressions which are functionally iden- 
tical to Eq. d?), but where (e.g. for temperature) the quantities 

Xj = T'i{r^f^, Tt = wri{wr),XT = ivr|2/(ivr|2) 

appear in place of Xb, Tb^Xb- In the following we will also use 
the unlabeled symbols X, T , x for, respectively, normalized 
fluctuations, fluxes and dissipation of a generic scalar. 

3.2. Scalar fluxes 

The normalized, conditional, averages of the temperature and 
salinity fluxes in our numerical experiments, E(Tt\Xt)/Xt and 
E(T^s \Xs)/Xs , are reported in Figs. |4la,b), computed in the pos- 
itive half-plane. The linear expression E{f\X)/X = X fits very 
well the data, particularly at larger Rayleigh numbers. Small 
deviations from linearity are evident only for large temperature 
or salinity fluctuations. Note that no further constants appear in 
fliis expression due to the identity (T) = / E{T\X)P(X)dX and 
the requirements {X^} = 1 and (T) = 1. This result has a very 
simple interpretation in the flow under study, if temperature and 
salinity are considered as almost passive scalars: as a parcel of 
fluid moves downward (upward) in the presence of background 
gradients, it carries with it the temperature and salinity values 
corresponding to a higher (lower) level, thus generating posi- 
tive (negative) temperature and salinity fluctuations. Molecular 
diffusion tends to remove the fluctuation, until an equilibrium 
between these two competitive factors is found, leading to a 
fluctuation proportional to the vertical speed of the parcel of 
fluid. From this, the relationship E{'F\X)/X - X trivially fol- 
lows. 

A diff'erent mechanism must apply for buoyancy fluctua- 
tions, which play a direct role in determining vertical accel- 
erations in the flow and which are generated by the doubly- 
difiiisive mechanism, rather than by vertical displacement of 
the fluid. Fig. lUc) reports the normalized conditional averages 
of buoyancy fluxes, E(Tb\Xb)IXb- While the linear expression 
E(Tb\Xb)IXb - Xb again fits the data, more severe deviations 
from linearity are evident at large buoyancy values, particularly 
at low Rayleigh numbers. This result is consistent with a sce- 
nario where the advective fluxes are determined to a large ex- 
tent by the motion of blobs of buoyancy of characteristic size 
I, travelling, on average, at a vertical velocity w determined by 
the balance between viscous drag and buoyancy forces. This 
idea is supported by the observation that the typical Reynolds 
numbers of the structures present in the buoyancy field are low 
and of order one (see table 1). If the blobs can be approximated 
as spheres of diameter /, carrying an average buoyancy fluc- 
tuation B' , the balance between the buoyancy and the Stokes 
drag forces would be RsB'f' oc Iw. If I is roughly indepen- 
dent of B' we may take B' oc w for uniform blobs. Note that 
the power- law exponents for the dependence on Rs of B' , w, 
Ix and in Table 1, are in good agreement with the relation- 
ship w oc B'fiRs , after taking the standard deviation of B' and 
w as the characteristic values for the buoyancy and the vertical 
velocity fields. In this simple scenario the advective buoyancy 
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Figure 4: (a) Plot of E{'Ft\^t)I^t vs. Xj for the four simulations at Rg = 
10^10',10'°,10"; data are gathered for one hundred convective times, once 
statistical stationarity has been reached. For clarity the latter three curves are 
shifted upward by one, two and three units, respectively. The straight lines are 
the theoretical fit E('F\X)/X = X. Only the positive normalized fluctuation 
half-plane is shown. For the negative half-plane odd symmetry holds, (b, c) 
The analogous plots for salinity and buoyancy. 
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Figure 5: Average vertical velocity vs. average buoyancy in connected regions 
with IB'I > IcTgi for the Rs = 10* and the Rs = lO" simulations. Errobars in 
gray are 95% bands obtained by jackknife subsampling. Data are gathered for 
one hundred convective times, once statistical stationarity has been reached. 



fluxes are proportional to the square of the buoyancy fluctu- 
ations, leading to E{Tb\^b)I^b - Xg- This scenario applies 
mainly at high Rayleigh numbers, where the organization of 
the flow in approximately spherical blobs is more marked. The 
internal structure of the blobs and their interactions will lead 
to deviations from linearity, which are apparent in Fig. Hfc), 
particularly at large buoyancy and at low Rayleigh number. To 
have a further insight into this issue we have partitioned the 
buoyancy field into connected regions with \B'\ > 2a b' and we 
have computed the average buoyancy and the average vertical 
velocity within each of these regions. The results are shown in 
Fig. |5] It is evident that at the lowest Rayleigh number the ver- 
tical velocity of these connected regions is fairly independent of 
the region's average buoyancy, thus breaking the proportional- 
ity w cc B' for extreme values of B' . At high Rayleigh number 
connected regions with a very high buoyancy move substan- 
tially faster than those having a lower average buoyancy. In this 
case the proportionality w oc B' , on average, holds fairly well. 

3.3. Scalar dissipation rates 

The conditional averages of temperature and salinity dis- 
sipation for the numerical experiments are reported in Figs. 
|6ja,b). To interpret them it is again useful to take temperature 
and salinity, to a good approximation, as passive scalars. The 
Kolmogorov-Obukhov-Corrsin scenario (e.g. (author?) 11211 ') 



assumes that the dissipation of variance of a passive scalar is, 
on average, independent of the concentration of the scalar it- 
self, which would yield a constant E(x\X) as a function of X. 
Sinai and Yakhot (author?) |11] suggested that concentration 
and dissipation of a passive scalar may be correlated. They 
modeled the dissipation as a parabolic function of concentration 
and linked this behaviour with the appearence of non-Gaussian 
tails. In our simulations the even parabolic expression 



fits well the temperature and salinity dissipation data (also in 
this case the normalization of X and ofx reduces the number of 
free constants in Eq. ([8])). Conditional dissipation is fairly con- 
stant at Rs - 10^ and smoothly assumes an upward parabolic 
shape as the Rayleigh number increases. At low Rs - 10^ we 
find a slightly downward shape: high fluctuations of tempera- 
ture or salinity are slightly less subject to dissipation than small 
ones. This may reflect the particular distribution of temperature 
and salinity inside blobs: extreme fluctuations are found mainly 
in the regions of low gradients at the core of the blobs, where 
they are protected from dissipation. 

Fig. |6tc) shows instead that at all Rayleigh numbers the dis- 
sipation of buoyancy is never independent of buoyancy itself. 
To some extent also the dissipation of buoyancy can still be in- 
terpreted as the contribution of a passive scalar. Any buoyancy 
fluctuation which is not aggregated into a blob is likely to be 
dynamically irrelevant: small buoyancy structures significantly 
different from rising or sinking blobs will be strongly damped 
by viscosity, and quickly dissipated. Therefore buoyancy in the 
background between the blobs is passively transported, and a 
parabolic shape of E(xb\Xb) at moderate values of Xb is to be 
expected. For large values of Xb the dissipation is dominated 
by the contribution of the blobs. The detailed distribution of 
buoyancy fluctuations within each blob and the distribution of 
blob sizes and intensities will determine the form of the tails of 
E{xb\Xb) (author?) d. At this stage we limit ourselves to 
report that the numerical simulations at high Rayleigh numbers 
show remarcably linear tails, and postpone to a future work an 
in-depth investigation of the morphology of the blobs. 

A simple fitting expression that allows us to join the picture 
inside and outside the blobs is 



E(xb\Xb) + 



aX\ 



(9) 



E(x\X)^\-y + yX'- 



(8) 



1 + b\XB\ 

The parameter k weights the Kolmogorean part of the dissipa- 
tion; a is the coeflicient of the quadratic component in the ex- 
pression of E{xb\Xb)\ ajb is the asymptotic slope of the linear 
tails. Although the constants k, a, b are functionally linked by a 
normalization constraint, writing this relationship explicitly is 
not as straightforward as in the case of the conditional expecta- 
tion of the fluxes. Here we prefer to independently fit all three 
constants that appear in Eq. (|9]l with a nonlinear least-square 
regression to the numerical data. The results are shown in Fig. 
|6jc) and show very good agreement. 

3.4. Agreement with the scalar distributions 

Since Eq. O and its equivalents for temperature and salin- 
ity represent exact expressions for the scalar fluctuation ampli- 
tude distributions, computing conditional averages of the scalar 
fluxes and dissipations from our experimental data and substi- 
tuting them into these equations, leads trivially to an almost 
exact overlap with the curves in Fig. |3](not shown). When a 
simple dependency E{T\X)IX - X \s used for the fluxes, to- 
gether with the fits ([8]l (for the dissipations of temperature and 
salinity) and (|9]i (for the dissipation of buoyancy), we obtain 
a very good agreement with the distributions computed from 
the numerical simulation data, as shown by the black lines in 
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Figure 6: (a) Conditional averaged dissipations of temperature and salinity, 
E(xt\Xt) and E(xs \Xs), for the four simulations at Rs = 10^ lO', 10'°, lO"; 
data are gathered for one hundred convective times, once statistical stationar- 
ity has been reached. For clarity the latter three curves are shifted upward by 
one, two and three units, respectively. Overimposed on the numerical data is 
the parabolical fit (8), where the constant y is determined by a nonlinear least- 
squares fit. Only the positive normahzed fluctuation half-plane is shown. For 
the negative half-plane even symmetry holds, (c) The analogous plot for buoy- 
ancy dissipation rates, with overimposed fit ^9). 



Figs. Oa-c). The only exceptions are the tails of the buoyancy 
distribution at low Rayleigh number, where the agreement is 
worse due to the deviations from linearity of E(Tb\Xb)I^b dis- 
cussed in Sec. 13.21 Note that, when E{T\X)IX -X'ls, assumed, 
the non-Gaussianity of the tails of the amplitude distributions 
is controlled by the form of the expected dissipation: the linear 
tails of Eq. ^ determine the exponential tails of the distribu- 
tions of buoyancy fluctuations, while a dissipation independent 
of scalar fluctuations (as shown by temperature and salinity at 
low Rs), leads to Gaussian distributions. 

4. Conclusions 

In the numerical experiments of fingering convection re- 
ported in this letter we find sharp evidence for exponential non- 
Gaussian tails in the buoyancy fluctuation distributions at high 
Rayleigh number In contrast, the statistics of temperature and 
salinity remain closer to Gaussianity even when those of buoy- 
ancy are already significantly non-Gaussian. As shown by us- 
ing a custom version of a theory by Yakhot (1989), this ob- 
servation can be understood in terms of the different properties 
of dissipation of a scalar directly creating vertical accelerations, 
such as buoyancy, compared to the dissipations of the individual 
buoyancy-changing scalars, such as temperature and salinity. 

There are some analogies with the phenomenology of 
Rayleigh-Benard convection, where, in the high Rayleigh num- 
ber "hard-turbulence" regime, temperature fluctuations present 
exponential-like tails. In that setting temperature fluctuations 
are equivalent to buoyancy, while in fingering convection this 
role is played by a linear combination of temperature and salin- 
ity. We suspect that the faster appearance, in fingering con- 
vection, of non-Gaussian statistics of buoyancy, compared with 
those of other active scalars, may apply also to other flows with 
multiple active scalars. 

The simple conceptual models introduced in Sees. (13.2b and 
(13.31 ) highlight the role of coherent dynamical structures in de- 
termining the vertical fluxes and dissipation of buoyancy fluc- 
tuations and consequently the appearance of non-Gaussian tails 
in their amplitude distribution. The same statistics for tempera- 
ture and salinity are closer to what could be expected of passive 
scalars. The blobs of fingering convection are generated by a 
very different mechanism compared to the plumes of Rayleigh- 
Benard convection and live in a small range of spatial scales 
where the effects of molecular diffusion and viscosity are very 
strong, but nonlinear terms are just as important. We have used 
a simple threshold in buoyancy to partition the flow and iden- 
tify the blobs but a better identification method will be needed 
to characterize in detail their structure and to explore their La- 
grangian properties, such as their mean free path, their charac- 
teristic life time and the dynamics of their mutual interactions. 

The changes in the shape of the distributions that we observe 
occur gradually over a range of Rayleigh numbers spanning 
three orders of magnitude. The convective fluxes and other in- 
dicators (Table [1]) follow cleanly simple scaling laws, with no 
evident breaks. This suggests that any change in the flow pat- 
terns aff'ecting the distributions is not as dynamically important 
as, for example, the changes in the plumes in the analogous 
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transition in Rayleigh-Benard convection. Nevertheless, we are 
quite confident that we are not yet observing any sort of ul- 
timate scaling regime of fingering convection. In fact, as we 
have argued in Sec. 13.21 the vertical convective fluxes can be 
modeled in terms of the equilibrium between a blob's buoy- 
ancy and a Stokes drag. However the scalings of Table[T]imply 
that the Reynolds number of an individual blob increases with 
the Rayleigh number, and, eventually, it will significantly ex- 
ceed one. At that point the dynamics will necessarily change, 
as the blobs will be subject to a drag having a nonlinear depen- 
dence on velocity. Whether this will simply mark a change in 
the slope of the scaling laws and in the form of the amplitude 
distributions, or if it will trigger more dramatic changes, such 
as the formation of the elusive staircases, remains to be seen. 
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